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We update our results for the heavy-heavy spectrum at /3 = 5.7, using NRQCD. This includes a scaling 
comparison with the T spectrum at /3 = 6.0 and first lattice predictions for the Be spectrum. 



1. Introduction 

We have recently given accurate results for bot- 
tomonium and charmonium spectroscopy on the 
lattice , exploiting the non-relativistic nature 
of these systems by using a non-relativistic effec- 
tive theory, NRQCD§. This can be systemati- 
cally matched to full QCD, order by order in 
and v'^/c^, where v is the velocity of the heavy 
quark in the bound state. A limited number of 
terms in the NRQCD action suffice provided that 
v'^ is small and the lattice heavy quark mass. Ma 
is not too small. The latter condition prevents a 
being taken to zero in this effective theory, but it 
should still be true that physical results scale in a 
limited region of a. This is all that is necessary to 
obtain continuum values. We test the scaling of 
our T spectrum results by comparing results at f3 
=5.7 with those at /3 = 6.0, an updated picture 
of which is presented here in the talks by Paul 
McCallum and Junko Shigemitsu. 

There is much current interest in the spec- 
trum of mixed bound states of bottom and charm 
quarks with the hope that they will soon be ob- 
served experimentally. Here we present results for 
Be spectroscopy from NRQCD and use our T and 
^ results to give good indications of the size and 
direction of systematic errors. We compare our 
results to recent calculations in potential models 
i 



2. Lattice NRQCD 

The NRQCD action, details of how heavy 
quark propagators are evaluated, the smearing 
functions and the fitting procedures used are de- 
scribed in detail in several publications [^[0J^ and 
are reviewed in this Proceedings in the talks by 
Paul McCallum and Junko Shigemitsu. 

The starting point of NRQCD is to expand the 
original QCD Lagrangian in powers of , the typ- 
ical quark velocity in a bound state. For the J/5' 
system ~ 0.3 and for T v"^ ^ 0.1. For Be 

~ 0.15 for the single particle with mass equal 
to the reduced mass of the 6, c system, but the 
kinetic energy will be shared unequally between 
the h and c. ~ 0.5 for the c quark so relativis- 
tic corrections will be more important than for cc 
systems. NRQCD enables these corrections to be 
added systematically. 

The action used here is the same one used 
in refs. |]l| and ||]. Relativistic corrections 
0{Mqv^) have been included for both h and c 
quarks. Other sources of systematic error are dis- 
cretisation errors and errors from the absence of 
virtual quark loops because we use quenched con- 
figurations generated with the standard plaquette 
action. Finite volume errors should be negligible 
because of the relatively small size of heavy-heavy 
systems. 

To calculate masses for be bound states we de- 
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fine b and c quark Green functions on the lattice, 
in the standard way. We used 200 quenched gluon 
field configurations at /3 = 6/(7^ = 5.7 generously 
suppfied by the UKQCD collaboration and 
fixed to Coulomb gauge, using a Fourier accel- 
erated steepest descents algorithm 0. The bare 
masses of the c and b quarks were chosen from fits 
to the kinetic mass of the 77c and the T, giving 
aruc = 0.8 (with stability parameter, n = A) and 
arrib = 3.15 (n = 2). 

We combine a 6 propagator with a c propagator 
to make states. From the evolution equations, 
it is more efficient to have numerous b quark prop- 
agators with small values of the stability parame- 
ter n and only one c quark propagator. We there- 
fore fixed the c quark with a 6 function source and 
used various smearing functions as sources for the 
b quark. This enabled us to also map out the T 
spectrum at the same time. 

Local meson operators are tabulated in [Q . Us- 
ing the notation we have looked at me- 
son propagators for the following states: ^5*0, '^Si, 
^Pi,^Po,^Pi and for both the E and T repre- 
sentation. Since C is not a good quantum number 
for the Be system the ^Pi and ^Pi mesons will 
mix and so in addition we calculated the cross- 
correlation between these two. For the s states, 
smearing functions both for the ground and first 
radially excited state were used as well as a local 
S function. From this all possible combinations 
of smearing at the source with a local sink were 
formed. For the p states only the ground state 
smearing function was used at the source. We 
calculated the dispersion relation for the T and 
for the (Be) by looking at the meson propagator 
for small non-zero momenta. Because of the rela- 
tively small size of these systems it is possible to 
use more than one starting site for the mesons. 
We used 8 different spatial origins and 2 different 
starting times to increase our statistics, but we 
bin correlation functions over each configuration. 

3. Simulation results 

The most accurate quantities we calculate are 
spin-independent splittings such as the spin- 
averaged IP-IS splitting. A comparison between 
this splitting at /? = 5.7 and (3=6.0 shows good 
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Figure 1. The scaling plot for the ratio nip/ (IP — 
lS)r at 3 different values of /?, filled circles. The 
nip values are preliminary UKQCD results using 
the tadpole-improved clover action. The T re- 
sults at large a are NRQCD results at /3 = 5.7 
and 6.0. The point at smallest a is taken from 
a UK(NR)QCD result using only the leading or- 
der NRQCD action. The open circle shows the 
experimental result. 



scaling behaviour once 0{a^) errors in the gluon 
field configurations are removed by estimating 
their effect on IS levels either perturbatively or 
non-perturbatively |^ (a 10% effect for the T at 
P =5.7). Figure 1 shows a scaling plot for the 
ratio of the IP- IS splitting to rUp. To see scaling 
we must use results for rUp which have 0{a) er- 
rors removed. The results in Figure 1 come from 
UKQCD |] using a tadpole-improved clover ac- 
tion for light propagators. The T splitting at /3 = 
6.2. is from ref. using a lowest order NRQCD 
action only and therefore has a potential 10% sys- 
tematic uncertainty. 

The difference with experiment for this ratio is 
similar to that we find when comparing splittings 
for bb and cc at (3=5.7 ^ and interpret it as a re- 
sult of incorrect running of the coupling constant 
between scales in the quenched approximation. 
The light hadron spectrum may still have large 
systematic errors arising from other sources such 
as finite volume so it is too early to put this result 
down entirely to quenching. 

We currently determine spin splittings only to 
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Figure 2. A comparison of T spin splittings for 
/3=5.7 and (3 — 6.0. The spin-average of s and p 
states is set to zero. 



leading order with no discretisation corrections. 
This has a large effect at f3 = 5.7 as can be 
seen from Figure 2. The discretisation corrections 
tend to reduce the overall size of splittings and 
also tend to equalise the Pi and ^Pi Po 

splittings. This is already evident aX 13 — 6.0 and 
is worse at j3 = 5.7. This can be improved by sys- 
tematically adding discretisation corrections to 
the E and B fields. 

Figure 3 shows results for the Be spectrum. We 
are able to predict the overall mass scale as well as 
the splittings since the c and h masses were fixed 
from the T and r^c masses. The fact that is 
different for cc and hh in the quenched approxima- 
tion means that we must fix for he separately. 
We do this by setting the IP-IS splitting to be the 
same as that for bb and cc (it is Mq independent). 
This gives an between the two previous val- 
ues and means that the quark masses in physical 
units are not the same in the three systems. This 
is an unavoidable artefact of quenching and adds 
to the systematic error. 

The Be mass is obtained by calculating correla- 
tion functions at non-zero momenta to derive the 
dispersion relation. The mass is the kinetic mass 
in the equation: 



E{p) = Po + pV2Mfe, 



(1) 



Mkin will not agree with the static mass, Po, be- 
cause the quark mass has been removed from the 
Hamiltonian. There is then a shift between Po 
and Mkin that must be applied per quark in a 
meson. We can check non-perturbatively that the 
shift for the Be is the average of that for cc and bb. 
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Figure 3. NRQCD simulation results for the spec- 
trum of the Be system using an inverse lattice 
spacing of 1.32 GeV. Error bars are shown where 
visible and only indicate statistical uncertainties. 



This is shown in Table 1 and works very well. For 
the clover action at large mass it does not work, 
and instead the shift is dependent, incorrectly, on 
the meson itself 

The comparison of our results with a potential 
model shows fairly good agreement (Figure 3) and 
indicates that relativistic effects may fortuitously 
not be very large. We have included the first rel- 
ativistic corrections but potential model calcula- 
tions include relativistic effects only through the 
phenomenological potential. Better estimates of 
these corrections need to be made. 

Our experience with heavy-heavy spectra in the 
quenched approximation indicates that our calcu- 
lation is likely to overestimate the 2S-1S splitting 
and to underestimate the hyperfine and p states 
splittings. Discretisation errors also reduce our 
spin splittings, as above. Unquenching and dis- 
cretisation corrections would then bring us into 
better agreement with the potential model results 
except for the p fine structure where we already 
have larger splittings. We believe as a result that 
the potential model calculation is underestimat- 
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Eo 


Mkin 


shift 


T 


0.5030 


7.04 


6.54 


Vc 


0.618 


2.430 


1.812 


Be 


0.6052 


4.79 


4.18 



Table 1 

Results for static and kinetic masses and the dif- 
ference between them for heavy-heavy mesons, 
using NRQCD at /3 =5.7. Values are in lattice 
units. 



ing these splittings. 

Our preliminary value for the Be decay con- 
stant is 460 MeV. This calculation includes only 
the leading order term in the heavy-light current. 
There is a 1/Mc correction for Be whereas in bb 
and cc the first correction appears at l/Mg. 

After this talk was given we learned that the 
ALEPH collaboration has candidate Be events 

il- 

4. Summary 

We present heavyonium spectrum results at f3 
= 5.7, including first lattice predictions for the Be 
spectrum. Future calculations will include further 
relativistic and discretisation corrections as well 
as a comparison to unquenched results at similar 
lattice spacing. 
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